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Abstract 

In SSVEP-based Brain-Computer Interface (BCI), it is very innportant to get an evoked EEG with a high signal to noise ratio 
(SNR). The SNR of SSVEP is fundamentally related to the characteristics of stimulus, such as its intensity and frequency, and it 
is also related to both the reference electrode and the active electrode. In the past, with SSVEP-based BCI, often the 
potential at 'Cz', the average potential at all electrodes or the average mastoid potential, were statically selected as the 
reference. In conjunction, a certain electrode in the occipital area was statically selected as the active electrode for all 
stimuli. This work proposed a dynamic selection method for the reference electrode, in which all electrodes can be looked 
upon as active electrodes, while an electrode which can result in the maximum sum relative-power of a specific frequency 
SSVEP can be confirmed dynamically and considered as the optimum reference electrode for that specific frequency 
stimulus. Comparing this dynamic selection method with previous methods, in which 'Cz', the average potential at all 
electrodes or the average mastoid potential were selected as the reference electrode, it is demonstrated that the SNR of 
SSVEP is improved significantly as is the accuracy of SSVEP detection. 
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Introduction 

SSVEP-based BCI system possesses many advantages compared 
to other types of BCI system [1,2,3,4,5,6], and one of the most 
prominent properties is its high transfer rate [7,8,9]. To get a high 
transfer rate, besides using a valid method for SSVEP extraction, it 
is most important to get an evoked EEG witli a high signal to noise 
ratio (SNR). In an SSVEP-based BCI system, a widely used 
method for SSVEP extraction is to compute the SSVEP relative- 
power by EFT [4,10,11], and this method is referred to Power 
Spectrum (PS) Method. In fact, the relative-power of SSVEP can 
be seen as the SNR of SSVEP in a specific frequency band. In the 
Power Spectrum Method, a reference electrode (for example, 'Cz') 
is frrsdy selected, and then one or a few active electrodes (for 
example, 'Ol' and/or '02') are selected [10,12]. The relative- 
power of a certain frequency in spontaneous EEG at the active 
electrode is computed to build a threshold, and the relative-power 
of the corresponding frequency in evoked EEG at the active 
electrode is computed within a short period, such as 2 s or 3 s, 
then compared with the threshold. If the relative-power of that 
frequency is higher than the corresponding threshold, it can be 
concluded that this frequency SSVEP is included in this span, in 
other words, it can be detected that the subject is staring at the 
button with this frequency flicker inside and decides to select that 
button. 

The BCI system which uses PS method is simple to build, 
because only one reference and one or a few active electrodes are 
employed [4,8,10,13,14,15,16]. However, because of the inter- 
subject dilference of the SSVEP power [17,18,19,20,21], a certain 



active electrode may be very effective at a known frequency 
stimulus for a few subjects, while not so valid for the other subjects 
[22,23,24,25]. Alternatively, even for the same subject, because of 
the traveling property of SSVEP [26], a certain active electrode 
may be very effective for a certain frequency stimulus while not so 
effective for the other frequency stimuli. In other words, the most 
effective active electrode can vary between subjects or different 
stimulus frequencies. To statically select only one or a few active 
electrodes, we cannot cover all of the most effective active 
electrodes, and thus must limit improvements to the SNR of 
SSVEP for all frequency stimuli. 

Consider inter-subject differences and the travehng property of 
SSVEP. If all the electrodes on the scalp are selected as the active 
electrodes and the relative-power of SSVEP at these electrodes are 
summed together as an SSVEP indicator, the drawback of only 
utilizing one or a few active electrodes can be overcome to some 
extent. Although different frequency SSVEPs can come to their 
maximum power at different electrodes, these maximum powers 
can all be included in the sum relative-power. The signal at each 
electrode can make contributions to the recognition of SSVEP 
frequency, so the detection accuracy can be improved for all 
stimuli. 

For a certain frequency SSVEP, the sum relative-power of 
SSVEP can vary with different reference electrodes [17]. In order 
to get the maximum sum relative-power, a suitable reference 
electrode should be confirmed dynamically for each frequency, 
and this reference electrode is referred to the optimum reference. 
The selection of the optimum reference electrode is conducted for 
every stimulus automatically. The optimum reference electrode 
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can vary between subjects and different stimulus frequencies. The 
sum relative-power under the optimum reference is used as an 
indicator of SSVEP. The method of dynamically selecting an 
optimum reference electrode under the situation of selecting all 
electrodes as active electrodes is proposed for the first time in this 
work, and is referred to as the Dynamic Selection (DS) Method. 

In this study, six frequencies in different bands were selected as 
the stimulus frequency, a 129 channel EEG system was used to 
record EEG signals, and 100 s length spontaneous EEG and 
evoked EEG were collected separately and then divided into 2 s 
length segments. For th(' e\'()k(;d EEG scgnw'nts, the sum relative- 
power of SSVEP for each frequency was computed according to 
the DS method, and compared to the sum relative-power under 
'Cz' reference, the reference of the average potential at all 
electrodes, or the reference of the average mastoid potential, 
respectively. The results indicate that the sum relative-power 
under the optimum reference is significantiy higher than the sum 
relative-power under other three kinds of reference; accordingly, 
the detection accuracy under the optimum reference is higher than 
that under the other three kinds of reference. Although the 
optimum reference electrodes for a certain frequency or a certain 
subject can be different from each other, most of them locate at the 
occipital lobe. 

There are many other kinds of BCI system except for SSVEP- 
based BCI, such as P300-based BCI [27], sensorimotor rhythm 
(SMR)-based BCI [28]. The singal used for classfication in BCI 
can mainly be recorded over the scalp, over the cortical surface, 
and within the brain [29]. A new classification method of BCI 
which is similar to the communication .system is proposed recentiy 
[29], according to this method, the BCI system can be sorted into 
five types: TDMA, FDMA, CDMA, SDMA, and HMA. SSVEP- 
based BCI is tiie kind of FDMA, while P300-based BCI is tiie kind 
of TDMA. In these different kinds of BCI system, although the 
signal extraction methods are different, for example, FFT method 
is often used in SSVEP-based BCI, while the superposition 
method is often used in P300-based BCI, it is very important to get 
an evoked EEG with a high SNR. Although the method proposed 
in this study is based on the property of SSVEP, it can be extended 
to other type BCIs. For example, in a motor imagery based BCI, 
the ERP amphtude is related to the reference electrode. To select 
dynamically an optimum reference can lead to a high ERP 
amplitude, which can improve the classfication accuracy of the 
BCI system. 

Methodology 

2.1 Ethics Statement 

This study was approved by the Human Research and Ethics 
Committee of the University of Electronic Science and Technol- 
ogy of China. Before the experiment, all the subjects were told the 
purpose and procedure of the experiment in detail and signed a 
consent form. These forms were approved by the University of 
Electronic Science and Technology of China Ethics Committee 
Data Acquisition. 

2.2 Data Acquirement 

Eleven subjects were chosen to take part in this experiment, 
having either normal sight or corrected normal sight. Six subjects 
were male, and five subjects were female. The mean age was 24 
(range 24±2) years old. The subjects were seated in a dark room, 
60 cm from the stimulator. A high luminance focused LED was 
used as the SSVEP stimulator, driven by a pulse generator, and 
the duty-cycle of the pulse was set to 1:1. The cycle of the pulse 
can only be adjusted in 1 ms step, and six cycles were selected, i.e. 



30, 40, 60, 80, 120, and 160 ms, the corresponding six frequencies 
were 33.33, 25, 16.67, 12.5, 8.33, and 6.25 Hz, located at P, a, 
and 9 band, respectively. This frequency arrangement can be used 
to study the validity of DS method for the fundamental frequency 
and the harmonics extraction. The order of these six stimuli was 
random for different subject, and this can cancel the probable 
influence resulted by the stimuli order. A 1 29-channel EEG system 
was used to collect the spontaneous EEG and evoked EEG. 
Electrode impedance was kept below 10 kfi, and salt water 
dropped into the electrode periodically in order to retain good 
contact with the subjects scalp. Fig. 1 shows the location of the 
electrodes in this system. In the EEG recording stage, 'Cz' 
electrode (No. 129) was selected as the reference. In order to avoid 
power line interference, the cutoff frequency of the EEG system 
was set to 49 Hz. For each stimulus, 100 s length spontaneous 
EEG was collected first, and then 100 s length evoked EEG was 
collected. This spontaneous EEG was used to buUd the threshold 
for each frequency SSVEP. 

2.3 Computation Method 

Because the SSVEP has relative immunity to noise such as eye 
or body movement [2,18], no pre-process method, such as 
removing eye movement was adopted in this work. 

2.3.1 SSVEP Gain Under Different Kinds of 
Reference. The original EEG was referenced at 'Cz' electrode. 
For 100 s length spontaneous EEG, FFT was applied directly, 
then a spectrum was attained at each electrode with a frequency 
resolution of 0.01 Hz. For a specific frequency 'f Hz and the 
electrode Li, L2, ... L„, the relative-power of 'f Hz at any 
electrode 'm' can be computed as follows: 

Rfm = Ptm/mean(P(f_i)^,P(f+i)„) (1) 

Where 'm' is from '1' to 'n', 'Rfm' stands for the relative-power of 
'f Hz at electrode 'm', 'Pfn,' stands for the absolute-power of 'f Hz 
at electrode 'm', 'mean(P(f !)„-,, P(f+i)ni)' stands for the average 
absolute-power from 'f-1' Hz to 'f+1' Hz at electrode 'm'. The 
sum relative-power of 'f Hz can be stated as: 

Rf = Rfi+RG-H...Rf„ (2) 

Where 'Rf' stands for the sum relative-power of 'f Hz at all 
electrodes in spontaneous EEG, this sum relative-power is looked 
as SNR of this frequency in spontaneous EEG, and it is used as the 
baseline of this frequency. In this work, 'n' equals to 1 29, and 'f 
equals to 33.33, 25, 16.67, 12.5, 8.33, and 6.25 Hz, respectively. 

For evoked EEG, the same method as that for spontaneous 
EEG was adopted, and the sum relative-power of SSVEP was 
obtained which is looked as SNR of this frequency in evoked EEG. 
This sum relative-power of SSVEP was divided by the 
corresponding baseline in spontaneous EEG and the quotient 
referred to as the SSVEP gain under 'Cz' reference: 

Gf = SRf/Rf (3) 

Where 'G,' stands for SSVEP gain at 'f Hz, 'SR,' stands for the 
sum relative-power of SSVEP at 'f Hz in the evoked EEG, i.e. the 
SNR of SSVEP in evoked EEG, and the computation of 'SRf' in 
evoked EEG is the same as that of 'Rf in spontaneous EEG. This 
SSVEP gain under 'Cz' reference is compared to the next 
computed SSVEP gain under the optimum reference. 

To compute the SSVEP gain under the other kinds of reference, 
the spontaneous EEG and evoked EEG at each electode were 
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Figure 1. Electrodes location of 129 channel EEG system. 

doi:10.1371/journal.pone.0104248.g001 



firstiy re-refrenced at the average potential at all electrodes and the 
average potential at both mastoids, respectively. The average 
potential at electrode No. 56 and 107 was viewed as the average 
mastoid potential. The computation of SSVEP gain under these 
two kinds of references is the same as that under the 'Cz' 
reference. 

2.3.2 Selecting the Optimum Reference. For 100 s length 
evoked EEG of a known frequency, firstly the channel No. 1 is 
selected as the reference and the original signal at each electrode 
re-computed under this reference to get a new evoked EEG. FFT 
is applied on the new evoked EEG. The sum relative-power of 
SSVEP is computed using the same method as in 2.3.1. Following 
this, channel No. 2 is selected as the reference, and the same 
method as above applied. These steps are repeated until aU 
channels have been selected as the reference once. 

For the 129 channel EEG system, a total of 129 sum relative- 
power is computed and the maximum within them chosen as the 
optimum sum relative-power of SSVEP. The corresponding 
electrode is then selected as the optimum reference for the known 
frequency. 

2.3.3 SSVEP Gain Under the Optimum Reference. For 

spontaneous EEG, the sum relative-power of each stimulus 
frequency is re-computed under the corresponding optimum 
reference, and viewed as the baseline of each frequency under the 
optimum reference. 

For evoked EEG, the sum relative-power of stimulus frequency 
is computed under the optimum reference, and this sum relative- 



power is divided by the corresponding basehne to get SSVEP gain 
under the optimum reference. The sum relative-power of the other 
five frequencies, which can be viewed as noise under the known 
stimulus, are computed under the corresponding optimum 
reference also, and divided by the corresponding baseline under 
the optimum reference to get noise gain. Noise gain is used to 
evaluate, while improving the SSVEP gain significantly under the 
optimum reference, whether the noise gain had improved 
appreciably or not. 

2.3.4 Detection Accuracy Under Different Kinds of 
Reference. Spontaneous EEG and evoked EEG were divided 
into segments of 2 s length separately. For each segment, '0' series 
were added with a length of 2 s to get a whole 4 s length segment. 
The technique of adding '0' series is utilized to improve the 
frequency resolution of FFT, and is widely applied in BCI studies 
[10,1 1,30]. In this work, if applying FFT direcdy on the 2 s length 
signal, the frequency resolution is 0.5 Hz. After appending '0' 
series, the frequency resolution can be improved to 0.25 Hz. Then 
the four kinds of reference, i.e. 'Cz' reference, the reference of the 
average potential at all electrodes, the reference of average 
potential at mastoid, and the optimum reference, are applied to 
this 4 s length segment respectively to compute the sum relative- 
power. The sum relative-powers extracted from each spontaneous 
EEG segment is used to build a threshold for each frequency. An 
estimated threshold is used to check the sum relative-power of a 
certain frequency in spontaneous EEG segments, and adjusted 
continually until the detection accuracy equaled 90%. This means 
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that the relative-power in 90% spontaneous EEG segments is 
smaller than this threshold and the adjusted value can be seen as 
the threshold of the corresponding frequency under a kind of 
reference. As a result, six thresholds can be confirmed under each 
type reference, respectively. 

Under different kinds of reference, the sum relative-power in 
evoked EEG segments is compared to the corresponding threshold 
to check whether the SSVEP is included. First, only the first 
harmonic is utilized for detecting SSVEP with the checking 
standard being as follows: For an evoked EEG segment including a 
known frequency SSVEP, if the sum relative-power of the first 
harmonic exceeded its threshold, while the sum relative-powers of 
other frequencies except the second harmonic are all below the 
corresponding thresholds, the detection for this segment is correct. 
Secondly, for the high frequency stimuli such as 33.33 and 25 Hz, 
only the first harmonic is used for detecting, while for other middle 
and low frequency stimuli, both the first and second harmonics are 
utilized for detecting SSVEP. The checking standard used was as 
follows: For an evoked EEG segment including a known frequency 
SSVEP, if the sum relative-power of the first harmonic or the 
second harmonic exceeded the corresponding threshold, while the 
sum relative-powers of other frequencies are all below the 
corresponding thresholds, then the detection for this segment is 
correct. 

For the results obtained via these steps, in order to test the 
significance of difference between methods, one-way Analysis of 
Variance (ANOVA) was applied. Significance level 'p' was selected 
as 0.05. If 'p' is smaller than 0.05, it suggests that there is a 
significant difference between the compared situations. 

Results 

3.1 Optimum Reference Distribution 

For the total 1 1 subjects, everyone took 6 SSVEP-frequency 
tests, so there were. 

66 optimum reference electrodes chosen and most of these 
electrodes were located at the occipital area. For a certain subject, 
the optimum reference for different frequencies can be different. 
For a certain stimulus frequency, the optimum reference for 
different subjects can vary. Table 1 shows the optimum reference 
for different subjects at different stimuli. Fig. 2 shows the optimum 
reference distribution topography for all subjects under all stimuli. 

3.2 SSVEP's SNR Improvement Under the Optimum 
Reference 

Under a 33.33 Hz stimulus, the SNR of SSVEP under the 
optimum reference is significantly higher than that under 'Cz' 
reference (F(l,20) = 4.45, p = 0.04), the average mastoid reference 
(F(l,20) = 6.8, p = 0.03) and the commom average reference 
(F(l,20) = 5.47, p = 0.03), while for the other 5 frequencies which 
can be looked at as noise, the ANOVA results ^p' are far bigger 
than 0.05. This suggests that there is no significant improvement of 
the SNR for noise under the optimum reference. The SSVEP's 
SNR improvement under the optimum reference can be seen from 
the average SSVEP gain across all subjects. For the stimulus 
frequency 33.33 Hz, the SSVEP gain under the optimum 
reference is 1.65 times of that under 'Cz' reference, 1.57 times 
of that under the average mastoid reference, and 1 . 7 times of that 
under the common average reference, while for other noise 
frequencies, this ratio is 1.0 or so, which means that the noise 
power in evoked EEG is the same level as that in spontaneous 
EEG under the optimum reference. 

Under a 25 Hz stimulus, the comparison results between the 
optimum reference and the other kinds of reference were similar to 
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that under 33.33 Hz stimuli. Under other middle or low frequency 
stimuli, i.e 16.67, 12.5, 8.33, and 6.25 Hz, the comparison results 
between the optimum reference and other three kinds of reference 
were similar to that under 33.33 Hz stimulus except that of the 
second harmonic. When stimulated at middle or low frequency, 
sometimes the second harmonic becomes stronger than the first 
harmonic under 'Cz' reference, or the average mastoid reference, 
or the common average reference, and the SNR of the second 
harmonic can be improved significantiy when utilizing the 
optimum reference. Table 2 lists the average sum relative-power 
(SNR) across all subjects and the ANOVA results under every 
stimulus. 

3.3 Detection Accuracy Under Different Kinds of 
Reference 

When only taking the first harmonic into account, the average 
detection accuracy across all subjects and stimuli is low under all 
kinds of references, i.e. 31.2% under 'Cz' reference, 34.4% under 
the average mastoid reference, 33.2% under the common average 
reference and 42.9% under the optimum reference. However, the 
average detection accuracy across all stimuli under the optimum 
reference is significantly higher than that under 'Cz' reference (F(l, 
20) = 8.19, p = 0.01), the average mastoid reference (F(l, 20) = 6.3, 
p = 0.01), and the common average reference (F(l, 20) =10.2, 
p = 0.01). For all subjects, the detection accuracy of different 
frequencies is different. Normally, the detection accuracy of 
middle and low frequencies such as 16.67, 12.5, 8.33, and 6.25 Hz 
is lower than that of high frequencies such as 33.33 and 25 Hz. 



Fig. 3 illustrates the average detection accuracy across all stimuli 
for every subject when only taking the first harmonic into account. 

When taking the first harmonic into account for high frequency 
stimuli such as 33.33 and 25 Hz (due to not recording the second 
harmonic), and taking the first and second harmonic into account 
for middle and low frequency stimuli, such as 16.67, 12.5, 8.33, 
and 6.25 Hz, the average detection accuracies across all stimuli 
can be improved significantly under all kinds of reference. 
Furthermore, the average detection accuracies across all subjects 
and stimuli are 62.2% under 'Cz' reference, 61.8% under the 
average mastoid, 61.6% under the common average and 73.7% 
under the optimum reference, respectively. In this situation, the 
average detection accuracy across all stimuli under the optimum 
reference is significantly higher than that under 'Cz' reference 
(F(l,20) = 1 1.3, p = 0.0), the average mastoid reference (F(l, 
20) = 9.52, p = 0.0), and the common average (F(l, 20) = 8.3, 
p = 0.0). For all subjects, the detection accuracy of different 
frequencies is different. Normally, the detection accuracy of 
middle and low frequency stimuli such as 16.67, 12.5, 8.33, and 
6.25 Hz is higher than that of high frequency stimuli such as 33.33 
and 25 Hz. Fig. 4 illustrates the average detection accuracy across 
all stimuli for every subject when taking both the first and second 
harmonic into account. 

Table 3 lists the average detection results across all subjects 
under different stimulus frequencies. From this table, it can be seen 
that when using middle or low frequency stimuli, the second 
harmonic is very important for SSVEP detection. Under any 
situation, the ANOVA 'p' is far smaller than 0.05, which suggests 
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Figure 3. Average detection accuracy across all stimuli for every subject when only taking the first harmonic into account. 

doi:10.1371/journal.pone.0104248.g003 
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Figure 4. Average detection accuracy across all stimuli for every subject when taking the first and second harmonic into account. 
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that the detection accuracy under optimum reference is signifi- 
candy higher than that under other three kinds of reference. 

Discussion 

In this work, even in the consideration of the second harmonic, 
the detection accuracy for some subjects is still low compared to 
other works [10,11], this is resulted by the experiment itself 
SSVEP power is related to factors such as stimulus intensity, 
frequency, modulation depth and to the subjects themselves. This 
work shows that there is considerable intra-difference for the same 
subject between different frequency stimuli, and there is great 
inter-difiference between subjects even under the same stimuli. 
These differences lead to the different detection accuracy ratios 
between subjects. The experiment in this work can lead to more 
inter-subject or intra-subject differences. Although every subject 
was asked to be seated in the same location, we have not 
confirmed whether the subject evoked a maximum SSVEP. 
Because a focused LED was used as the stimulator of SSVEP in 
this study, if the subject did not stare at the flicker from the correct 
vision angle, the light projecting into the eyes decreased acutely 
and the maximum SSVEP was not evoked. In fact, we have 
studied the spectrum of the 100 s length evoked EEG, in some 
trials, the peak is not clear, which suggests that the SSVEP is very 
weak. Therefore, detection accuracy in this work is not as high as 
that in the other real BCI experiments [10,1 1]. However, this does 
not matter for the comparison between the methods for reference 
electrode selection. In a real SSVEP-based BCI, the experiment 
design is very important. Firstiy, suitable frequency arrangements 
can lead to a high SSVEP power and reduce the interference 
between stimuli. For example, normally selecting a middle or low 
frequency can evoke a large SSVEP, and a certain stimulus 
frequency should be far enough from other stimuli and their 
harmonics. Secondly, stimulator selection is also very important. 
We have compared the influence on SSVEP power by different 
stimulators [31]. When using a CRT monitor or emanative LED, 
the subject can receive almost the same strength light in a wide 
range, so a litde shift of the vision angle does not affect SSVEP 
power. However, when using a focused LED as the stimulator, the 
light concentrates in a narrow area, the shift of the vision angle has 
a great influence on SSVEP power. So, in a multi-target BCI 
system, because of the vision angle shifting significantly, the 
focused LED should be avoided. 

In order to check the influence on the second harmonic by the 
dynamic selection method, the SNR of the second harmonic was 
also computed Firstly, only the first harmonic was used to detect 
SSVEP. Then the first and second harmonic were used together to 
detect SSVEP. The results show that, when stimulating at a 
middle or low frequency, the SNR of the second harmonic can be 
improved significantly when using the optimum reference, 
accordingly the detection accuracy using the first and second 
harmonic together is significandy higher than when only using the 
first harmonic. In order to avoid the power line interference of 
.50 Hz, the cutoff frequency of the recording system is set 49 Hz, 
so a high frequency harmonic of stimulus such as 25 and 33.33 Hz 
was not collected. In consequence, the detection accuracy for these 
stimuli is smaller than that of the low frequency stimuli. In an 
SSVEP-based BCI, the harmonics are very important for 
improving detection accuracy, so the first and second harmonics, 
sometimes even the third harmonic, should be taken into account 
even under the optimum reference. In fact, some works have 
discussed the usefulness of harmonics in detail [4, 1 0] . 

In past SSVEP-based BCI systems, often one or a few active 
electrodes were chosen from the electrodes located at the occipital 
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area [1,10,1 1,13,32]. There are some drawbacks to doing this. For 
the same frequency, the optimum active electrode in the occipital 
area may be different because of inter-subject differences. In order 

to get the highest detection accuracy for each subject, it should be 
tested first to confirm the optimum active electrode before 
detecting targets. This however, is a waste of manpower and time 
[11,24,33]. Although the occipital area is considered the source of 
SSVEP [25], suggesting SSVEPs in this area are normally higher 
than those in other areas, for different frequencies, the location of 
the maximum SSVEP in the occipital area is different. In a multi- 
targets SSVEP-based BCI system, there ar(; normally many 
frequencies applied, it is impossible to find an electrode in the 
occipital area which is optimum for all frequencies. So, ail the 
electrodes are selected as the active electrode in this work. Under 
this selection, for any subject and any stimulus frequency, there is 
no need to confirm the optimum active electrode before detection. 
The SSVEP power in different areas can all be collected as an 
indicator of SSVEP, and this indicator includes more information 
than that for only one or a few active electrodes selected. 
Therefore, the detection accuracy for all stimuli is improved 
significandy. 

When selecting all electrodes as the active electrode, it is 
important to find an optimum reference under which the sum 
relative-power of SSVEP is at a maximum. Under the optimum 
reference, the SNR of SSVEP is improved significantly compared 
to the other kinds of reference. While using the optimum reference 
to improve the SSVEP gain, if the noise gain is improved the same 
level as SSVEP gain, optimum reference makes no improvement 
for SSVEP detection accuracy. Fortunately, except for the second 
harmonic, the other noise relative-power is not increased under 
optimum reference. Consequentially, when using the first and 
second harmonic together to detect SSVEP under the optimum 
reference, a higher detection accuracy compared to the other 
kinds of reference can be attained. 

Although the optimum references for different subjects or 
different stimulus frequencies can differ, most of them are located 
at the occipital area. This can be understood through the following 
analysis. The SSVEP is a response to the visual stimulus mainly by 
the primary visual cortex, and the occipital area is the source of 
SSVEP. Here we hypothesize an ideal situation, i.e. there is an 
ideal reference, which is not related to EEG at any electrode, and 
only one electrode 'Eo' at the occipital area is the SSVEP source. 
Furthermore, the SSVEP at electrodes 'Ei',... 'Ej,,'... 'E^' are 
transferred from 'Eq' via the scalp. The SSVEP between the ideal 
reference and electrode 'Eq', 'Ej',... 'E^'... 'E^' are expressed: 



these amplitudes, the initial ®o> ®i> ••• ®m-.. ®n can be different 
also. 

When using electrode 'Eq' as the reference, the relative SSVEP 
at other electrodes 'Ei',..- 'Em'... 'En' can be stood by: 

Si =Ai * sm{(o * t-l-0i) — Ao * sin((B * t+(j>o) 



Sm = Am * sin(co * t-l-i^„)-Ao * sm{(o * t+ (j)^) (5) 



Sn = An * sin(a) * t-l-(^n) * sin((B * t+(j)o 



And the sum power Tq' of these signals is: 



Po = ^* ^ ([M * sin(a» * t-l-i/ii) — Aq * sin(a) * t-l-i^o)]^ + - 



-I- [An * sin(a) * t-l-(^n) "-^o * sin(<» * i + <l>o)] )dt 



(6) 



Where 'T' is the cycle of the frequency 'cd'. 

When using electrode 'E,,-,' as the reference, the relative SSVEP 
at other electrodes is represented: 



50 = Ao * sin(a) * t-l-(^o) — A„ * sin(a) * t + (l)„ 

51 = Ai * sin(a) * t-l-(^i) — A„ * sin(a) * t + (l)„ 

Sn=An * sin(ft) * t+<j)^) — A„ * sin((B * t + (j)„ 
And the sum power 'P^' of these signals is. 



([Ao * sin(a) * t + ^!io)-A„ * sin(co * t + </>„)] 



(7) 



(8) 



-I- [An * sin(co * t+(j)„) —A„ * sin(co * t+<l>„)] )dt 



The 'Pq' and 'Pm' can be computed, and the difference of the 
sum power under reference 'Ep' and 'E^' is illustrated as: 



(4) 



Ao * sin(co * t + (^o) 
Ai * sin((o * t-l-0i) 

Am * sin(a) * t-l-(^„) 

An * sin((B *t+<j>^) 



A(), Ai,.An,...A„ are the SSVEP amplitude at different 
electrodes compared to the ideal reference, 'ra' is the SSVEP 
frequency, <l>o, ^i, ...<!>„,■.■ 'tn ^ue the initial phase of SSVEP at 
different electrodes. Because of the travelling property of SSVEP 
and the attenuation characteristics of the scalp, the amplitude Aq, 
Ai, ... Am, ... An can be different and Aq is the biggest among 



Po—Pm = 



(9) 



{n-l)*Al (n-\)*Al 

2 2 
—A\ * Ao* cos((^i — <^o) 

— A2 * Aq* COS(<j>2—<j>o) — -- 

-A„*Ao*cosi(l>„-(l>o) 

+ Ao* A„*cos((l>o-<t>„,) 
+ Ai *A„* cos(0, + . 
+A„*A„*cos((l>„-(t>J 



This difference is related to SSVEP amplitude and initial phase 
at every electrode. The different attenuation properties of a 
subject's scalp can lead to a different distribution of SSVEP. Even 
for the same subject, the distribution of SSVEP amplitude and 
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initial phase can be different for different stimulus frequencies. It is 
impossible to prove 'Po-Pm' is bigger or smaller than zero using 
mathematical methods. In order to understand this, a group data 
including 129 amplitudes and 129 initial phases are simulated 
randomly 100,000 times in a computer to compute equation (9). 
Except for assuming the source electrode 'Eo' with the maximum 
amplitude, there is no other limitation for other amplitudes and all 
initial phases, this suggests these simulant data can include all the 
real SSVEP amphtude and initial phase. The results show that the 
difference 'Po-P,-,,' is bigger than zero in more than an average 
probability of 98'%. In other words, even using the real amplitudes 
and initial phases in place of those in equation (9), 'Pq' is bigger 
than 'P„i' in most situations. This suggests that, if using SSVEP 
source as the reference, a maximum sum relative-power can be 
attained. Therefore, in the dynamic selection method, most 
optimum references locate at the occipital area. 

There is a very great difference between the method proposed 
in this work and the method in which th(- a\ (;rage potential at all 
electrodes in the time domain are selected as the reference. When 
selecting the average potential at all electrodes in the time domain 
as a reference, although the background noise can be canceled to 
some extent, more amplitude information of SSVEP has been 
canceled on average because of the traveling features of SSVEP, 
and thus cannot lead to the maximum sum relative-power of 
SSVEP. In this work, the sum rclati\T-power at all electrodes in 
the frequency domain is selected as the indicator of SSVEP. The 
phase of SSVEP has not been taken into account, and the SSVEP 
amplitude information at all electrodes remains in the sum 
relative-power. Maximum sum relative-power of SSVEP can be 
identified. 

The cost of the dynamic selection method is the increased 
complexity of the system because more electrodes can be lixed. 
Except for this drawback, the procedures using the dynamic 
selection method are ea.sy to apply. During the stage of building a 
threshold for each stimulus, the spontaneous EEG within a time 
period is first obtained, then each stimulus frequency is used once 
to evoke SSVEP for a period of length the same as for the, 
spontaneous EEG. Ultimately, the optimum reference for each 
stimulus is confirmed automatically, and the threshold under this 
optimum reference is computed automatically also. Therefore, 
time consumption in this stage is minimal. Alternatively, with the 
other rcft'rc'nce selection method, normally finding a suitable 
active electrode for all stimuli and confirming the threshold is very 
time consuming. In the formal detection stage, the procedure is 
similar to that of the power spectrum method, i.e. the power of 
every adopted frequency under the corresponding optimum 
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reference is computed and compared to the corresponding 
threshold. If every power is smaller than its corresponding 
threshold, then there is no button selected. If there is only one 

power bigger than its corresponding threshold, the button with the 
corresponding frequency flicker inside is selected. If there are more 
powers bigger than their corresponding thresholds, other tech- 
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Conclusion 

Compared to other SSVEP extraction methods, in which one 
reference is selected statically and one or a few electrodes in the 
occipital area are chosen as the active electrode. Dynamic 
Selection Method uses more active electrodes, and thus can 
increase the complexity of the system. The method of Dynamic 
Selection Method improves SSVEP's SNR and detection accuracy 
significantiy and is easy to employ by decreasing the number of 
electrodes, thus being applicable in a real time SSVEP-based BCI. 
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